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ABSTRACT 

The  Jordan  sorting  problem  is,  given  the  intersection  points  of  a  Jordan 
cur\e  with  the  x-axis  in  the  order  in  which  they  occur  along  the  curve,  to  sort 
them  into  the  order  in  which  they  occur  along  the  x-axis.  This  problem  arises 
in  clipping  a  simple  polygon  against  a  rectangle  (a  ‘window")  and  in  efficient 
algorithms  for  triangulating  a  simple  polygon.  Hoffmam  Mehlhom.  Rosen- 
stiehl,  and  Tarjan  proposed  an  algorithm  that  solves  the  Jordan  sorting  problem 
in  time  linear  in  the  number  of  intersection  points,  but  their  algorithm  requires 
the  use  of  a  sophisticated  data  structure,  the  level-linked  search  tree.  We  pro¬ 
pose  a  variant  of  the  algorithm  of  Hoffman  et  al.  that  retains  the  linear  time 
bound  but  simplifies  both  the  operations  required  on  the  key  data  structure  and 
the  data  strocture  itself. 
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1.  Polygon  Clipping  via  Jordan  Sorting 

Let  F  be  a  simple  n-vertex  planar  polygon  and  let  W  be  a  rectangular  window  in  the 
plane.  The  problem  of  clipping  P  against  W  is  to  compute  the  polygon  or  set  of  polygons  S 
that  bound  the  intersection  of  the  interior  of  P  with  the  interior  of  W.  The  polygon  clipping 
problem  is  fundamental  in  computational  geometry  and  in  computer  graphics;  it  arises,  for 
example,  in  window  management  and  in  hidden  surface  removal  [12].  In  the  laner  application, 
the  problem  is  actually  more  general,  in  that  both  P  and  W  can  be  polygons  with  holes,  but  our 
results  on  the  basic  problem  extend  to  the  more  general  problem.  In  our  discussion,  we 
assume  that  every  intersection  point  of  P  and  W  is  a  crossing  point,  i.e.,  a  point  where  P 
crosses  from  the  inside  to  the  outside  of  W.  Our  approach  extends  easily  to  handle  portions  of 
P  that  are  tangent  to  W  in  points  or  line  segments 

The  process  of  polygon  clipping  divides  naturally  into  two  steps: 

(1)  Find  the  points  of  intersection  of  P  and  W. 

(2)  Find  the  polygonal  curves  into  which  the  points  of  inteisection  divide  P  and  W,  and 

group  these  curves  to  form  the  set  S  of  output  polygons. 

Step  1  involves  primarily  geometric  computation.  It  can  be  carried  out  in  O(n)  time  by 
proceeding  from  vertex  to  vertex  along  P  and  determining,  for  each  line  segment  of  P,  whether 
and  where  it  crosses  W.  An  asymptotic  running-time  bound  of  O(n)  follows  immediately  from 
the  fact  that  there  are  O(n)  intersection  points.  Minimizing  the  constant  factor  in  the  running 
time  involves  interesting  issues;  see,  e.g.,  [7].  This  method  of  performing  Step  1  produces  the 
intersection  points  in  the  order  in  which  they  occur  aiong  P,  and  as  a  side  effect  determines  the 
polygonal  curves  into  which  the  intersection  points  divide  P. 

Step  2,  on  the  other  hand,  involves  primarily  topological  computation.  The  hard  pan  of 
the  computation  is  to  son  the  intersection  points  into  the  order  in  which  they  occur  alone  W. 
Once  this  is  done,  producing  the  set  of  output  polygons  is  straightforward  [12]. 

Much  work  on  polygon  clipping  has  concentrated  on  Step  1  at  the  expense  of  Step  2. 
Sutherland  and  Hodgman  [9]  dealt  with  the  sorting  problem  correctly  (in  an  appendix),  but  pro¬ 
posed  a  grouping  method  that  is  based  on  the  computation  of  shortest  paths  and  that  need  not 
produce  the  correct  output  polygons.  Liang  and  Barsky  [5]  ignored  the  need  for  Step  2 
entirely.  Weiler  and  Atherton  [12]  described  a  correct  grouping  method  but  did  not  explicitly 
discuss  the  sorting  problem. 

If  a  general  sorting  algorithm  is  used  u.»  sort  the  intersection  points,  the  time  for  Step  2. 
and  hence  for  the  entire  polygon  clipping  problem,  is  0(n  log/i).  From  an  algorithmic  point  of 
view,  the  interesting  aspect  of  polygon  clipping  is  that  a  general  sorting  algorithm  is  not 
required.  Indeed,  Step  2  can  be  performed  in  O(n)  lime,  and  hence  so  can  polygon  clipping 
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Let  us  abstract  the  sorting  problem  slightly  to  make  it  easier  to  study.  Consider  a  Jordan 
curve  C  in  the  plane  that  crosses  the  x-axis  everywhere  it  touches  it.  (A  closed  Jordan  curve  is 
a  homeomorphic  image  of  a  circle.  An  open  Jordan  curve  is  a  homeomorphic  image  of  a  line 
segment.  Most  topologists  use  the  term  “Jordan  curve”  to  mean  a  closed  Jordan  curve.  We 
extend  the  term,  however,  to  mean  either  an  open  or  a  closed  Jordan  curve;  the  key  property  of 
Jordan  curves  is  that  they  do  not  have  self-crossings.  To  be  concrete,  one  can  think  of  C  as  a 
simple  polygonal  curve.)  Let  zt,  z2,  •  •  •  .  zn  be  the  sequence  of  intersection  points  of  C  with 
the  x-axis  in  the  order  in  which  they  occur  along  C.  The  Jordan  sorting  problem  is  to  sort  the 
given  sequence  z\,  z2,  .  .  .  ,  z„  by  x-coordinate. 

An  algorithm  for  Jordan  sorting  can  be  applied  directly  to  the  sorting  part  of  Step  2  of 
polygon  clipping  Jordan  sorting  is  also  needed  in  efficient  algorithms  for  triangulating  a  sim¬ 
ple  polygon  [1,111.  Hoffman,  Mehlhom,  Rosenstiehl,  and  Taijan  [4]  devised  an  0(n)-time  Jor¬ 
dan  sorting  algorithm.  The  heart  of  their  method  is  a  sophisticated  data  structure,  the  level- 
linked  search  tree. 

Our  purpose  in  this  paper  is  to  improve  the  algorithm  of  Hoffman  et  al.  by  simplifying  it 
while  preserving  the  O(n)  time  bound.  First,  we  simplify  the  kinds  of  operations  needed  on 
the  key  data  structure.  We  show  that  Jordan  sorting  can  be  reduced  to  performing  an  inter¬ 
mixed  sequence  of  two  kinds  of  operations  on  a  collection  of  sorted  lists: 

(i)  insert  an  item  at  the  front  or  the  back  of  a  given  list; 

(ii)  find  an  item  in  a  given  list  and  split  the  list  just  before  or  just  after  this  item. 

The  algorithm  of  Hoffman  et  al.  requires  a  third  list  operation  in  which  a  middle  section 
of  a  list  is  split  out  and  the  two  end  sections  are  concatenated.  Since  our  algorithm  does  not 
require  this  third  list  operation,  we  obtain  a  second  simplification:  the  use  of  heterogeneous 
finger  trees  to  represent  the  lists,  in  place  of  level-linked  trees  (which  implement  homogeneous 
finger  trees). 

The  remainder  of  the  paper  consists  of  three  sections.  In  Section  2,  we  reduce  the  Jordan 
sorting  problem  to  the  list  manipulation  problem  discussed  above.  In  Section  3  we  show  that 
the  list  manipulation  problem  can  be  solved  in  O(n)  time  using  heterogeneous  finger  trees  or 
any  equivalently  efficient  data  structure.  Section  4  contains  some  concluding  remarks. 

2.  Jordan  Sorting  via  List  Manipulation 

Let  us  review  the  approach  of  Hoffman  et  al.  [4]  to  the  Jordan  sorting  problem.  (See  also 
[11].)  For  notationai  concision,  we  interpret  the  expression  “z,  <  z/'  to  mean  that  the  abscissa 
of  z,  is  less  than  the  abscissa  of  z;;  that  is,  we  consider  die  points  z  to  be  ordered  by  their  x- 
coordinates.  We  write  <z(_),z1>  to  denote  the  part  of  C  that  goes  from  z,_j  to  z,  without 
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crossing  the  x-axis;  thus,  if  i  h  j  mod  2,  <z;_i,z;>  and  <z;_i,z;>  lie  on  the  same  side  of  the  pr¬ 
axis.  We  also  write  <z(. .  .  .  ,zy>  to  denote  {j.^<<zk,zk+i>. 

The  sequence  zj,  z2 . z„  gives  rise  to  two  forests,  as  follows.  Assume  that  <zi,z2> 

lies  above  the  x-axis.  (If  not,  reflect  C  about  the  r-axis.)  For  1  <  /  <  n  let  /,  =  min  {z,_i,z.} 
and  r,  =  max  {z.^.z,}.  We  say  that  a  pair  {z,_i,z,}  encloses  a  point  z  if  /,  <  z  <  rt.  We  say 
that  a  pair  {z.-.j.z,}  encloses  a  pair  {z;_i,zy}  if  i  =  j  mod  2  and  {z,_i,z;}  encloses  both  z;_i 
and  zy;  thus,  <z,_i.zy>  lies  between  <z,_i,z,>  and  the  x-axis.  For  any  two  pairs  { z,_ i , zt }  and 
{z7_i,z;}  such  that  i  =  j  mod  2,  the  simplicity  of  C  implies  that  each  pair  encloses  an  even 
number  of  points  (zero  or  two)  of  the  other.  The  Hasse  diagram  of  the  “encloses”  relation  on 
the  set  of  pairs  { { z  2; - 1  >*2,}  I  1  ^  i  ^  Ln  / 2J)  is  a  forest,  called  the  upper  fore  ■  Similarly,  the 
Hasse  diagram  of  the  “encloses”  relation  on  the  set  of  pairs  {{*£>*21  +  1}  I 
1  Si  ^  L(/t-l)/2j}  is  called  the  lower  forest.  In  both  of  these  forests,  we  order  each  set  of 
siblings  by  placing  {z>_j ,z,}  before  {z,_),zy}  if  tt  <  ly  this  makes  each  forest  into  an  ordered 
forest  By  adding  a  dummy  pair  {-oc.oo}  to  each  forest,  we  create  two  ordered  trees,  called  the 
upper  tree  and  the  lower  tree.  (See  Figure  1.) 

A  family  is  a  set  of  pairs  consisting  of  a  pair  and  all  of  its  children  (in  the  appropriate 

tree).  Suppose  that  the  set  of  Finite  pairs  in  the  family  is  {z(i_i  ,z(| },  (zla_i  ,z,a } . 

{■Vi  -  A,  )•  If  the  family  does  not  include  the  pair  {-00,00},  then  it  corresponds  to  a  simple 

closed  curve  formed  from  <z;l_i,rll>,  <z,1_i,z,1> . <Al-i>z„>  along  with  appropriate  parts 

of  the  x-axis  interconnecting  them;  the  region  inside  this  curve  is  the  family  region.  If  the  fam¬ 
ily  includes  the  pair  {-00,00},  then  it  corresponds  to  a  simple  open  curve  formed  from 

<Zil_i,z,1>,  <Zij_i,zlj> . along  with  parts  of  the  x-axis  interconnecting  them  and 

two  rays  that  go  from  the  leftmost  curve  to  -00  and  from  the  rightmost  curve  to  00;  the  con¬ 
nected  region  bounded  by  this  curve  that  contains  no  other  parts  of  the  x-axis  is  the  family 
region.  The  family  region  for  a  family  in  the  upper  tree  lies  above  the  x-axis,  while  the  family 
region  for  a  family  in  the  lower  tree  lies  below  the  x-axis.  (See  Figure  2.) 

The  Jordan  sorting  algorithm  proceeds  incrementally,  processing  the  points  zj,  z2 . 

z„  one  at  a  time  and  building  the  upper  tree,  the  lower  tree,  and  a  list  of  the  points  in  sorted 
order.  Processing  the  point  z,  involves  inserting  z,  into  the  appropriate  position  in  the  sorted 
list  and,  if  1  >  1,  adding  {z,_i,z,}  to  the  appropriate  tree  (the  upper  tree  if  i  is  even,  the  lower 
tree  if  i  is  odd.) 

Before  examining  the  details  of  this  algorithm,  we  make  a  new  observation  about  the 
family  regions  that  simplifies  the  required  processing.  Suppose  that  z,,  z2,  .  .  .  ,  z,  have  been 
processed  so  far,  i.e.,  the  upper  and  lower  trees  have  been  constructed  representing  the  part  of 
curve  C  ir.-m  z,  to  z,.  Lei  F  be  the  family  so  far  constructed  containing  the  pair  {z,_]  ,z,}  as  a 
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child,  and  iet  R  be  the  corresponding  family  region;  observe  that  no  part  of  <zx . z,>  can 

lie  in  the  interior  of  R.  Let  z}  *■  z,_i  be  the  point  among  the  points  in  pairs  of  F  such  that  z, 
and  zj  are  adjacent,  i.e.,  no  point  in  a  pair  of  F  lies  between  z,  and  zy.  Consider  the  point  ?k 
where  C  first  enters  R  after  passing  through  z,  (if  C  ever  enters  R).  Then  either 

(a)  zk  lies  between  z,  and  zy  or 

(b)  the  parent  pair  in  F  is  {-»,«>}  and  z,  and  zk  bracket  ail  finite  points  in  pairs  of  F  other 

than  z,. 

(See  Figure  3.) 

We  prove  this  locality  lemma  by  contradiction.  Suppose  it  were  not  true.  Then  there 
would  be  two  finite  points  in  pairs  of  F,  say  zp  and  zq,  such  that  zp  but  not  zq  lies  between  z, 
and  zk.  (See  Figure  4.)  Since  R  is  a  path-connected  region,  there  is  a  simple  curve  S  lying 

entirely  inside  R  connecting  z,  and  z*.  Curve  S,  together  with  <z, . z*>,  forms  a  simple 

closed  curve  O  with  zp  on  its  inside  and  zq  on  its  outside.  The  part  of  <zj . zk>  that  con¬ 

nects  zp  and  zq  must  intersect  O.  Since  R  is  a  family  region,  this  intersection  cannot  occur  on 
S,  so  it  must  occur  on  <z„  .  .  .  ,z*>;  but  this  violates  the  simplicity  of  C.  This  contradiction 
establishes  the  lemma. 

The  locality  lemma  allows  us  to  represent  each  family  by  one  ot  two  lists,  called  sibling 
lists.  Consider  a  family  F  whose  most  recently  added  pair  is  {z,_i,z(}.  If  {z,_i,z,}  is  the 
parent  pair  of  F,  then  F  is  represented  by  one  sibling  list  containing  all  the  pairs  in  F  except 
{z,_i,z,}.  If  (z/_| ,Zj>  is  not  the  parent  pair  of  F,  then  F  is  represented  by  two  sibling  lists,  one 
containing  all  pairs  in  F  both  of  whose  points  are  less  than  or  equal  to  z,  and  the  other  contain¬ 
ing  all  pairs  in  F  both  of  whose  points  are  greater  than  or  equal  to  z,  (one  of  these  lists  may  be 
empty).  Sibling  lists  never  include  parent  pairs;  each  finite  pair  occurs  in  two  families  (once  as 
a  parent  and  once  as  a  child)  but  in  only  one  sibling  list. 

Now  we  are  ready  to  describe  the  details  of  the  Jordan  sorting  algorithm.  Initialization 
consists  of  creating  a  sorted  list  containing  ->»,  zlt  z2,  z 3,  and  »  in  sorted  order,  and  creating 
two  singleton  sibling  lists  containing  {zi,z2}  and  {z2,zj}.  The  main  part  of  the  algorithm 
consists  of  repeating  the  following  steps  for  i  from  4  to  n: 

Step  1.  Let  v  be  the  point  in  the  sorted  list  preceding  z,_j.  If  i  is  odd  and  v  =  zj,  replace  v  by 
its  predecessor  in  the  sorted  list.  (Point  z\  is  in  no  pair  in  the  lower  family.)  Let  {zy_|,zy}  be 
the  pair  containing  v  such  that  i  m  j  mod  2.  [Thus,  <zy_ltzy>  and  <z,_i,z,>  lie  on  the  same 
side  of  the  x-axis.]  (If  v  =  let  {zy_i,zy}  be  {-oo,oo}.) 

Step  2.  Let  w  be  the  point  in  the  sorted  list  following  z,_i.  If  i  is  odd  and  w  =  z\,  replace  w 
by  its  successor  in  the  sorted  list.  Let  (z*_j  .z* }  be  the  pair  containing  w  such  that 
i  *  k  mod  2.  [Thus,  <z*_i,z*>  and  <z,_i,z,>  lie  on  the  same  side  of  the  x-axis.]  (If  w  = 
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let  {2*_!,z*}  be  {-oo.oo}.) 

Step  3.  Assume  that  z,_,  <  z,;  the  other  case  is  symmetric. 

(a)  [Insert  {z,-_ltZj}  into  its  proper  place  with  respect  to  its  siblings  to  the  left.]  If  {z^j.z,} 
encloses  z,_j,  create  a  new  singleton  sibling  list  containing  {zJ_3  ,z, };  otherwise,  insert 
{z;_t,z,}  after  {z;_ltz;}  in  the  sibim.  list  containing  {z^_1(zy}.  (The  locality  lemma 
implies  that  {z^_i,z;}  is  last  on  this  list.) 

(b)  [Split  off  any  children  of  {z,_|,z,}.]  If  {z*_|,z*}  does  not  enclose  z,-_j,  split  the  sibling 
list  containing  {z*_i,z*}  into  two  lists,  one  containing  all  pairs  whose  points  are  less  than 
z;,  the  other  containing  all  pairs  whose  points  are  greater  than  z,.  (The  locality  lemma 
implies  that  {z*_i,z*}  is  first  on  its  sibling  list  before  the  list  is  split.)  One  or  the  other 
of  the  split  lists  may  be  empty. 

(c)  [Insert  z,  into  its  proper  place  in  the  sorted  list.]  If  {z,_i,z,}  has  no  children,  insert  z, 
immediately  after  z,_i  in  the  sorted  list.  Otherwise,  let  {zm_;,zm}  be  the  rightmost  child 
of  {z,_i,z,};  insert  z,  immediately  after  rm  in  the  sorted  list 

(See  Figure  5.) 

The  correctness  of  this  algorithm  follows  from  the  locality  lemma  and  other  properties  of 
Jordan  curves.  Let  us  make  a  few  comments  about  the  data  structures  needed  to  implement  the 
algorithm.  We  assume  the  input  is  given  in  the  form  of  a  doubly  linked  list,  with  next(z)  and 
prev(z)  being  the  successor  and  predecessor  of  intersection  point  z  along  C,  respectively.  The 
output  is  also  represented  by  a  doubly  linked  list,  with  after{z )  and  before(z)  being  the  succes¬ 
sor  and  predecessor  of  intersection  point  z  along  the  x-axis.  In  addition,  a  bit  even(zi )  indicates 
for  each  intersection  point  z,  whether  i  is  even  or  not.  With  this  representation,  it  is  not  neces¬ 
sary  to  construct  the  pairs  explicitly.  A  pair  {z,_1(z(}  can  be  represented  just  by  z,_ j,  since 
z,  =  next(zi- 1)  is  computable  in  constant  time.  Each  of  the  sibling  lists  is  then  just  a  list  of 
intersection  points,  rather  than  a  list  of  pairs  of  points. 

The  total  space  requirements  of  the  algorithm  are  four  pointers  and  one  bit  per  point,  plus 
whatever  space  is  necessary  to  represent  the  sibling  lists.  The  running  time  ot  the  algorithm  is 
O(n)  plus  the  time  needed  to  perform  O(rt)  insertion  and  split  operations  on  sibling  lists.  In 
the  next  section  we  discuss  a  data  structure  for  representing  the  sibling  lists  that  requires  O(n) 
time  and  space  for  the  list  operations. 

The  Jordan  sorting  algorithm  can  be  extended  easily  to  two  related  problems:  the  recogni¬ 
tion  problem,  i.e.,  testing  whether  a  given  sequence  can  arise  as  the  sequence  of  intersection 
points  of  a  Jordan  curve  with  the  x-axis  (see  [4]),  and  the  problem  of  Jordan  sorting  with 
error-correction  required  by  the  efficient  triangulalion  algorithm  of  Tatjan  and  Van  Wyk  [11]. 
Neither  extension  affects  the  0(n )  time  bound. 


-  7  - 


3.  Insertion  and  Splitting  Operations  on  Sorted  Lists 

We  have  reduced  the  Jordan  sorting  problem  to  a  pure  problem  in  data  structures,  which 
we  can  formulate  as  follows.  Maintain  a  collection  of  sorted  lists  whose  items  are  selected 
from  a  totally  ordered  universe  U,  subject  to  the  following  three  kinds  of  operations: 

(i)  make-lisKx ):  Create  a  new  singleton  list  containing  x. 

(ii)  instrt{x,y ):  If  x  <  y,  insert  x  at  the  front  of  the  list  containing  y;  in  this  case,  y  must 
initially  be  at  the  front  of  its  list.  If  x  >  y,  insert  x  at  the  back  of  the  list  containing  y, 
in  this  case,  y  must  initially  be  at  the  back  of  its  list. 

(iii)  split{x,L ):  Split  list  L  into  two  lists,  one  containing  all  items  less  than  or  equal  to  x 
and  the  other  containing  all  items  greater  than  x. 

A  sequence  of  m  intermixed  operations  is  to  be  performed  on  an  initially  empty  collection 
of  lists.  We  seek  a  list  representation  such  that  the  total  time  for  m  operations  is  0(m). 

A  suitable  list  representation  is  the  heterogeneous  finger  tree.  This  is  a  balanced  search 
tree  in  which  the  pointers  along  the  two  ribs  of  the  tree  (the  paths  from  the  root  to  the  first  and 
last  nodes  in  symmetric  order)  go  up  instead  of  down.  We  omit  a  detailed  description  of  the 
data  structure,  since  it  can  be  found  elsewhere  [11].  Heterogeneous  finger  trees  support  the 
desired  list  operations  in  the  following  amortized  time  bounds:  0(1)  for  make-list  and  insert, 
and  0(log(min  {l-k,k}  +  2))  for  split,  where  /  is  the  size  of  the  input  list  and  k  and  l-k  are 
the  sizes  of  the  two  output  lists.  By  “amortized”  time  bounds  we  mean  that  for  any  sequence 
of  operations,  the  sum  of  the  amortized  time  bounds  is  an  upper  bound  on  the  sum  of  the 
actual  times.  See  Taiwan’s  survey  paper  [10]  for  a  thorough  discussion  of  amortization. 

The  amortized  time  bound  for  split  is  too  large  to  conclude  immediately  that  the  total 
time  for  a  sequence  of  m  operations  starting  with  no  lists  is  0(m).  By  using  an  extra  amortiza¬ 
tion  argument,  however,  we  can  reduce  the  amortized  time  of  split  to  0(1).  The  reasoning  is 
the  same  as  that  used  by  Goldberg  and  Taijan  [3]  to  charge  splitting  time  to  concatenations. 
(In  our  case  this  time  is  charged  to  insertions.) 

We  use  the  idea  of  a  potential  function  [10].  We  define  the  potential  of  a  list  of  size  1  to 
be  c  (/  -  log  /),  where  c  is  the  constant  in  the  amortized  time  bound  for  splitting,  and  the  base 
of  the  logarithm  is  two.  We  define  the  total  potential  of  a  collection  of  lists  to  be  the  sum  of 
their  potentials,  and  the  nominal  time  of  a  list  operation  to  be  its  amortized  time  bound  plus 
the  net  increase  in  potential  it  causes.  For  any  sequence  of  list  operations,  the  sum  of  the  nom¬ 
inal  times  equals  the  sum  of  the  amortized  time  bounds  plus  the  final  total  potential  minus  the 
initial  total  potential.  The  initial  total  potential  is  zero  (there  are  no  lists  initially)  and  the  final 
total  potential  is  nonr.egative.  Thus  the  total  time  required  to  perform  any  sequence  of  m  list 
operations  is  at  most  the  sum  of  the  nominal  time  bounds. 
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The  nominal  time  to  initialize  a  list  is  0(1).  The  nominal  time  for  an  insertion  is  also 
0(1),  because  its  amortized  time  is  O  (1)  and  the  increase  in  potential  when  the  list  grows  from 
length  /  to  l  +1  is 

c(l  +1  -  log(/+l)  -  /  +  log/) 

=  C(1  +  log-1-) 

<c  =  0(1). 

Finally,  the  nominal  time  for  a  split  is  0(1);  when  sublists  of  sizes  k  t  1  and  l-k  >  1  are 
formed  by  a  split,  the  amortized  time  is  c(log  (min  { k,l-k }  +  2)),  and  the  potential  increase  is 

c(k  -  log  k  +  (l-k)  -  log  (l-k)  -1  +  log  /) 

=  c(-  log*  -  log (/— *)  +  log/) 

=  c(—  iogmin  {*,/-*}  +  log - rf  ---- ) 

max  {k,l-k} 

<.  c (-log min  {k,l-k}  +  1), 

since  /  <  2  max  {k,l-k}.  Therefore  the  sum  of  the  amortized  time  and  the  potential  increase  is 
at  most  3c  =  O  (1).  Since  the  nominal  time  for  each  kind  of  list  operation  is  0(1),  the  total 
time  for  m  list  operations  is  0(m),  as  desired. 

4.  Remarks 

There  are  two  obvious  ways  to  apply  the  Jordan  sorting  algorithm  to  clip  a  polygon 
against  a  convex  window.  Sutherland  and  Hodgman  [9]  suggest  that  the  polygon  be  clipped 
sequentially  against  each  line  that  contains  a  side  of  the  window;  the  set  of  polygons  resulting 
from  clipping  against  each  side  becomes  the  input  for  clipping  against  the  next  side.  It  may, 
however,  be  more  efficient  to  clip  the  polygon  against  the  entire  window  at  once.  To  do  this, 
we  need  only  be  able  to  compare  two  intersection  points  with  respect  to  their  order  around  the 
window  boundary.  If  we  choose  a  point  on  the  window  boundary  as  the  origin,  a  point  on  the 
interior  as  the  center,  and  a  direction  in  which  the  order  increases  (say  clockwise),  then  we  can 
perform  each  comparison  in  0(1)  time  by  computing  the  angles  of  the  rays  joining  the  center 
to  the  points  to  be  compared. 

Throughout  this  paper,  and  particularly  in  the  preceding  paragraph,  we  have  ignored  the 
time  required  to  find  the  points  of  intersection  between  the  n-sided  polygon  P  and  the  window 
W.  If  W  has  a  fixed  number  of  sides,  then  all  of  the  intersection  points  can  be  found  in  O  (n) 
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time,  by  a  generalization  of  the  method  described  in  Section  1.  If  W  is  convex  and  has  k  sides, 
then  we  can  traverse  P  to  find  the  intersection  points,  using  a  binary  search  algorithm  on  the 
sides  of  W  to  discover  whether  each  vertex  lies  inside  or  outside  the  polygon  and  when;  each 
edge  of  P  intersects  V',  the  total  time  for  this  method  is  O  (n\ogk).  If  W  is  not  convex,  then 
the  number  of  intersections  /  between  P  and  W  is  0  ( nkf,  the  algorithm  of  Chazelle  and 
Edelsbrunner  [1]  can  be  used  to  find  all  /  intersections  in  0  (/  +  (n  +k) log  (n  +k))  time,  and  in 
fact  produces  enough  information  to  perform  Step  2  directly. 

To  summarize,  when  the  clipping  window  W  has  a  fixed  number  of  sides,  Step  1  can  be 
performed  in  O  (n)  time,  so  the  complexity  of  Step  2  de'ermines  the  complexity  of  any 
polygon  clipping  algorithm.  When  W  has  a  variable  number  of  sides  or  is  not  convex,  how¬ 
ever,  the  dominance  of  Step  2  in  the  complexity  of  a  polygon  clipping  algorithm  is  no  longer 
so  obvious. 

As  noted  in  the  paper  by  Hoffman  et  al.  [4],  the  use  of  the  splay  tree  (8),  a  self-adjusting 
form  of  search  tree,  gives  a  very  simple  algorithm  for  Jordan  sorting  that  may  run  in  linear 
time.  This  method  may  well  be  preferable  in  practice  to  the  algorithm  we  have  presented,  even 
though  the  linear  time  bound  for  the  method  based  on  splay  trees  remains  only  a  conjecture. 
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Figure  Captions 


Figure  1.  A  Jordan  curve  ar. .  its  upper  and  lower  family  trees.  In  this  and  the  following  three  figures, 
the  intersection  z,  of  the  curve  with  the  x-axis  is  labelled  i. 

Figure  2.  The  family  region  corresponding  to  {24,25}  and  its  children  in  the  lower  family  tree  of  the 
curve  in  Figure  1. 

Figure  3.  The  locality  lemma  determines  where  C  can  re-enter  region  R  after  it  leaves  at  2,.  At  the  top, 
region  R  is  finite,  and  C  must  re-enter  it  between  z,  and  zr  At  the  bottom,  region  R  is  infinite,  and  C 
may  re-enter  it  to  the  right  of  z,  or  to  the  left  of  ad  finite  points. 

Figure  4.  In  this  illustration  of  the  proof  of  the  locality  lemma,  the  boundary  of  region  0  is  dashed. 
The  portion  of  the  boundary  outside  R  is  <2,,2)+1>, ....  <2t_1,2t>. 

Figure  5.  Illustration  of  Step  3  of  the  algorithm.  In  (i)  and  (ii),  substep  (a)  creates  the  singleton  sublist 
{2,_i,r,};  in  (iii)  and  (iv),  substep  (a)  adds  {2,-1, 2,}  after  {2y_t,2;}.  In  (i),  (lii),  and  (iv),  substep  (t) 
splits  off  zero  or  more  children  of  {r, ,2, }. 


